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Abstract — Codes for rank modulation have been recently 
proposed as a means of protecting flash memory devices from 
errors. We study basic coding theoretic problems for such codes, 
representing them as subsets of the set of permutations of n 
elements equipped with the Kendall tau distance. We derive 
several lower and upper bounds on the size of codes. These 
bounds enable us to establish the exact scaling of the size of 
optimal codes for large values of n. We also show the existence 
of codes whose size is within a constant factor of the sphere 
packing bound for any fixed number of errors. 

Index terms — Bose-Chowla theorem, flash memory, inversion, 
KendaU tau distance, rank permutation codes. 

I. Introduction 

Codes in permutations form a classical subject of coding 
theory. Various metric functions on the symmetric group (S„ 
have been considered, giving rise to diverse combinatorial 
problems. The most frequently studied metric on (3„ is the 
Hamming distance. Codes in (3„ with the Hamming distance, 
traditionally called permutation arrays, have been a subject of 
a large number of papers; see, e.g., the works by Blake et al. 
[1] and Colbourn et al. [5]. 

In this paper we are interested in a different metric on 6„ 
which we proceed to define. Let a = {<j{1), ... ,a{n)) be 
a permutation of the set [n] — {1,2, ... ,n}. The Kendall 
tau distance dr {<j,tt) from a to another permutation vr is 
defined as the minimum number of transpositions of pairwise 
adjacent elements required to change a into tt. Denote by 
Xn — (Snjdr) the metric space of permutations on n 
elements equipped with the distance dr. 

The Kendall distance originates in statistics and has been 
adopted as a measure of quality of codes under the so- 
called rank modulation scheme first considered by Chadwick 
and Kurz [3]. In this scheme, the transmitted sequences are 
given by permutations of n elements while information is 
carried by the relative magnitude (rank) of elements in the 
permutation rather than by the absolute value of the elements. 
The motivation for considering this scheme in [3] stems from 
systems in which transmitted signals are subjected to impulse 
noise that changes the value of the signal substantially but 
has less effect on the relative magnitude of the neighboring 

Alexander Barg is with tlie Department of Electrical and Computer Engi- 
neering and Institute for Systems Research!, University of Maryland, College 
Park, MD 20742 and Institute for Problems of Information Transmission, 
Moscow, Russia (e-mail: abarg@umd.edu). 

Arya Mazumdar is with the Department of Electrical and Computer 
Engineering and Institute for Systems Research, University of Maryland, 
College Park, MD 20742 (e-mail arya@umd.edu). 

Research supported in part by NSF grants CCF0830699, CCF0635271, 
DMS0807411. 



signals. Recently (and independently of [3]) rank modulation 
was suggested by Jiang et al. [14], [15] as a means of efficient 
writing of information into flash memories. Rewriting the 
contents of a group of memory cells is easy if one needs 
to increase the charges of the cells or leave some of them 
unchanged and impractical if some of the charges need to 
be decreased. Furthermore, reliability of the data stored in 
flash memory is affected by the drift in the charge of the cells 
caused for instance by aging devices or other reasons. Since 
the drift in different cells may occur at different speed, errors 
introduced in the data are adequately accounted for by tracking 
the relative value of adjacent cells, i.e., the Kendall distance 
between the groups of cells in memory. These considerations 
make rank modulation suitable for coding for flash memories. 
More details of both the writing and the error processes in 
memory are given in [14] and references in that paper 

The focus of our work is on bounds and constructions 
of codes in the Kendall space X„. Coding-theoretic con- 
siderations call for estimating the volume of the sphere in 
Xn because it can be used to derive basic bounds on the 
size of codes. Spheres in the Kendall space were studied by 
analytic means in a number of earlier works [17], [18] relying 
on the well-known correspondence of permutations and their 
inversion vectors; however it turned out that code bounds that 
can be obtained from these works do not lead to nontrivial 
(other than or 1) values of the code rate. Regarding specific 
code families for correcting Kendall errors, the only previous 
work is that by Jiang et al. [14] who constructed a family of 
single-error-correcting codes of size M > ^{n — 1)!, i.e., at 
least half the maximum possible. 

Our results. In this paper we discuss several possible ways 
to bound the size of codes for rank modulation of a given 
distance, often calling them rank permutation codes. We derive 
a Singleton-type bound and sphere -packing bounds on such 
codes. Since the maximum value of the distance of the distance 
in Xn is (2), this leaves a number of possibilities for the 
scaling rate of the distance for asymptotic analysis, ranging 
from d — 0{n) to d — Q{n^). These turn out to be the 
two extremes for the size of optimal rank permutation codes. 
Namely, earlier work in combinatorics of permutations implies 
that a code with distance d = O(n^) occupies a vanishing 
proportion of the space Xn while a code of distance 0{n) 
can take a close-to-one proportion of its volume. We cover the 
intermediate cases, showing that the size of optimal codes with 
distance d ~ < e < 1 scales as exp(nlnn(l — e)). It is 
interesting that unlike many other asymptotic coding problems, 
the Kendall space of permutations affords an exact answer for 
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the growth rate of the size of optimal codes. The proof of the 
bounds rehes on weight-preserving embeddings of Xn into 
other metric spaces which provide insights into the asymptotic 
size of codes. 

We also show the existence of a family of rank permutation 
codes that correct a constant number of errors and have size 
within a constant factor of the sphere packing bound. The 
construction relies on the well-known Bose-Chowla Theorem 
in additive number theory. 

Section II of our paper is devoted to the relation of 
the Kendall metric space to other metric spaces related to 
permutations. In Section III we use these insights to derive 
bounds on codes for rank modulation, and conduct their 
asymptotic analysis. Section IV contains a construction of t- 
error-correcting rank permutation codes. 

II. Weight-preserving embeddings of the Kendall 

METRIC SPACE 

We begin with recalling basic properties of the distance 
dr such as its relation to the number of inversions in the 
permutation, and weight-preserving embeddings of 6„ into 
other metric spaces. Their proofs and a detailed discussion 
are found for instance in the books by Comtet [6] or Knuth 
[16, Sect. 5.1.1]. 

The distance dr is a right-invariant metric which means 
that dr{(Ti,a2) = dr{(Tia, <72(t) for any ct, cti, (T2 G &n where 
the operation is the usual multiplication of permutations. This 
implies a number of simple properties: 

(a) we can define the weight of the permutation a as its 
distance to the identity permutation e = (1, 2, . . . , rt). 

(b) dr{-K,a) = d^(7r-\cr-l). 

Because of the invariance, the graph whose vertices are 
indexed by the permutations and edges connect permutations 
one Kendall step apart, is regular of degree n — 1. At the 
same time it is not distance-regular, and so the machinery of 
algebraic combinatorics does not apply to the analysis of code 
structure. The diameter of the space Xn equals N — (2) and 
is realized by pairs of opposite permutations such as (1,2,3,4) 
and (4,3,2,1). 

The main tool to study properties of dr is provided by 
the inversion vector of the permutation. An inversion in a 
permutation cr e 6„ is a pair {a{i), a{j)) such that i < j and 
> cr(j). It is easy to see that dr{(T,e) = ^(c), the total 
number of inversions in cr. Therefore, for any two permutations 
a-i,(T2 we have dr(o'i,o'2) = I{a2(Ji^) — ^i^i'^'i'^)- I" other 
words, 

dr{a,TT) = G [nY : i 7^i,7r(i) > TT{j),a{i) < a{j)}\. 

To a permutation a E (3„ we associate an inversion vector 



e G„ = Z2 X 



Z„, where x^{i) = |{j : j < 



i + I, > a{i + 1)}|, i — 1, . . . , n — 1 and is the 
set of integers modulo m. It is well known that the mapping 
from permutations to the space of inversion vectors is one-to- 
one, and any permutation can be easily reconstructed from its 
inversion vector Moreover, 



For the type of errors that we consider below we introduce 
the following £1 distance function on G„ : 

ri-l 

d(a;,y) = ^|x(^)-2/(^)| (a^,yeG„) (2) 

1=1 

where the computations are performed over the integers, and 
write 1 1 a; 1 1 for the corresponding weight function (this is not 
a properly defined norm because G„ is not a linear spacefl 
For instance, let <ti — (2, 1, 4, 3), 1T2 — (2, 3, 4, 1), then Xa-^ — 
lOljXcrj = 003. To compute the distance dr{ai,(J2) we find 

1(^2^1-1) = /((1, 4, 3, 2)) = 11 (0,1, 2)11 =3. 

Observe that the mapping cr ^ 0;^ is a weight-preserving 
bijection between X„ and the set G„. At the same time, since 
the groups S„ and G„ are not isomorphic (one is commutative 
while the other is not), this mapping is not distance-preserving. 
However, a weaker property is true, namely. 



dr(cri,Cr2) > d{Xcri,Xa 



(3) 



Indeed, transposing two neighboring entries of a permutation 
cr changes the inversion count I{a) by one, so the mapping 
Xn Gn preserves distances to the identity permutation. 
Thus, if there exists a code in G„ with £1 distance d then 
there exists a code in Xn with Kendall distance at least d. 

Another embedding of Xn is given by mapping each permu- 
tation to a binary iV-dimensional vector a whose coordinates 
are indexed by the pairs C [n]^, i < j, and a(i,j) = 1 if 
the pair {i,j) is an inversion and a(, = otherwise. Clearly 
the Hamming weight of a equals /(cr), and so this mapping is 
an isometry between Xn and a subset of the Hamming space. 
This mapping was first considered in [4]. 

III. Bounds on the size of rank permutation codes 

An (n, M, d) code C C Xn is a set of M permutations in 
which any two distinct permutations are at least d distance 
units apart. Let A{n, d) be the maximum size of the code in 
Xn with distance d. For the purposes of asymptotic analysis 
we define the rate of a code C C X„ of size M as R{C) ~ 
i'^*^ Let 

\nA{n,d) 



Inn! 



lim 

n — ^00 



In n\ 



(1) 



4=1 



be the capacity of rank permutation codes of distance d (as 
usual, we assume that the limit exists and use liminf and 
lim sup as appropriate if it does not). The main result of this 
section is given in the following theorem whose proof is given 
in Sections UlLBl and UlLCl below. 
Theorem 3.1: 

{1 if d = 0{n) 

1-e if fi = e(ni+"), < e < 1 (4) 
ifd = e(n2). 

Remark. As will be seen from the proof, the equality '^{d) — 
1 — e holds under a slightly weaker condition, namely, d — 
n^+^a(n), where a{n) grows slower than any positive power 
of n. 

' This metric is reminiscent of the asymmetric distance function on the set 
of integer n-tuples, p{x,y) = Yl"Zi(^{^) - [^l- 
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A. A Singleton bound 

Theorem 3.2: Let d > n ~ 1. then 



A{n, d) < [3/2 + \/n{n - 1) - 2d + i/ij !. 



(5) 



Proof: Let C be an {n, M, d) code. Since the metric dr 
is right invariant, we can assume that C contains the identity 
permutation e. 

Let k < n and G 6^ be a code derived from C in the 
following way. Let (/jj. : S„ 6/j be a mapping that acts on 
a by deleting elements k + 1, . . . ,n from it. Thus, ^^(cr) is a 
permutation on k elements that maintains the relative positions 
of the elements of [k] given by a. 

Let k be the greatest number such that (pk is not injective. 
Then (pk+i is injective, and M < {k + 1)1. Suppose that 
permutations <Ti,a2 S &n are such that ^^(cti) = 4'k{o'2)- 
Because of the last equality, none of the k symbols in (72 erf ^ 
contain pairs that form inversions. Therefore, 

fn\ fk^ 

d < dr(ai,C72) < 



This gives 



k < 



1 + v/4n(n- 1) - 8d+ 1 



which proves inequality (|5]). This estimate is nontrivial if 
I + y/ n{n — 1) — 2d + 1/4 < n which is equivalent to the 
condition d > n — 1. ■ 
To gain an insight into this bound, let d = SN. Using the 
inequality m! < (to/2)™ in (|5]l, we obtain the asymptotic 
inequality 



A{n, d) < exp(n 



As we will show in the next section, the — 5 in this bound 
can in fact be improved to o(l). 

B. Sphere packing bounds 

Denote by B,. = Br{Xn) the ball of radius r in the Kendall 
space Xn. Clearly, 



nl , , , n! 



(6) 



The embeddings of X„ into other metric spaces considered 
in the previous section can be used to derive estimates of 
A{n, d) based on these inequalities. In particular, computing 
the volume of the metric ball in G„ and using ([3]), we will 
derive a lower bound in (|6]l. At the same time, both lower and 
upper bounds will follow from the embedding of Xn in the 
Hamming space described above. 

Let Kn{k) — \{a ^ S„ : /„((t) = fc}| be the number of 
permutations with k inversions. By Kn{k) is the number 
of solutions of the equation 



n-l 

E 

i=l 



Xi = k, where Xi S Z.^+i. 



Then clearly Kn{k) — for k > N and 

Kn{k) = Kn[N -k^ for < fc < ]^N. 



The number of inversions in a random permutation is asymp- 
totically Gaussian with mean ^A^ and variance ^" w 
[10, p. 257]. The generating function for the numbers 
Kk{n) has the form 



K{z) 



k=0 



Kn{k)z' 



n 



(7) 



For 1 < k < n the number of permutations with k inversions 
can be found explicitly [16]: 



Knik) = 

J>2 



k 
k 
hfc 



u 

k — u 



n + A: — 3 
k-2 
1\ fn 



J-1 



k - Uj - j 



(8) 



where Uj = — j)/2 and the summation extends for as 
long as the binomial coefficients are positive (it contains about 
l.GVk terms). 

For 1 < k < n the expression for Kn{k) is given above. In 
particular, it implies that = n, and 3) < (n — 1)!. 
As shown in [18], for n = fc + to, to — > 00, k>Q 



Kn{k) = (0.289. 



om+n — 1 

.) ^ (1 + 0(to-1)). (9) 



The case of fc > n is much more difficult to analyze. An 
obvious route for finding asymptotic approximation of Kn{k) 
is to start with the integral representation of the coefficients 
of K{z) ©. Namely, since K{z) converges for every z in the 
finite plane, we can write 



Kn{k) 



1 

27ri 



n 

£=1 



1 - z 



dz. 



where C is a circle around the origin. Asymptotic analysis of 
this expression involves saddle point calculations and is rather 
involved [17]. The next theorem is a combination of results of 
Margolius [18] and Louchard and Prodinger [17], stated here 
in the form suitable for our context. 

Theorem 3.3: There exist constants ci and C2 such that 



Kn{k) < exp(cin) 
Knik) ~ n\/ exp(c2n) 



if fc = 0(n), 
if fc = Q{n^). 



The implicit constants in this theorem can be found in cited 
references. 

From this theorem and inequalities (|6]l, we obtain the two 
boundary cases of the expression for ^^{d) in Q. 

C. Bounds from embedding in the ii space 

In this section we prove the remaining case of Theorem l3.1l 
Our idea is to derive bounds on "^(d) by relating the Kendall 
metric to the £1 metric on (5„. From the results of Diaconis 
and Graham [9], 



V2-D ((71,(72) < dr(«7l,(72) < L'(ai,(72) 



(10) 



where D{ai,a2) = ^"=1 |o'i(*)~f2(«)|- Therefore, any code 
C C &n with Kendall distance d must have ti distance at least 
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d and any code C C 6„ with £i distance d must have Kendall 
distance at least d/2. 

Remark. Define T{ai, 02) to be the number of inversions of 
(not necessarily adjacent) symbols needed to change a\ into 
(72 • Paper [9] in fact shows that 

dr{cri,a2) + T{ai,a2) < D{ai,a2) 

which is a stronger inequality than the one given above. We 
however will not use it in the derivations below. 

Proposition 3.4: Let Br{Hn, x) be the metric ball of radius 
r with center at x in the space H„ = {l,2,...,n}" with 
the £1 metric. Then the maximum size of a code in Xn with 
distance d satisfies 
i! 

< A{n,d) < 



To prove the lower bound, write Bj.{z) as z + Dr{z), where 
Dr{z) is the set of differences: 

n 

Driz) = {m e Z" : < n - 1, 1 < i < n; ^ < r 

i=l 

and z + M e Hn}. 
Writing Br{l) in the same way as 1 + Z)+, we have 

n 

£»+ = {m e Z" : < < n - 1; ^ \u^\ < r}. 



maxccej?„ \B2d-iiH„,x)\ 



where t = [(d- 1)/2J. 

Proof: Under the trivial embedding S„ 7J„ the £1 
distance does not change, so any code C in 6„ with £1 distance 
d is also a code in Hn with the same distance and as such, 
must satisfy the Hamming bound. Together with (fTOl i this gives 
the upper bound of our statement. 

Turning to the lower bound, let us perform the standard 
"Gilbert procedure" in the space of permutations with respect 
to the £1 distance, aiming for a code V with £1 distance m. 
The resulting code satisfies 

\V\ max \B„,-iie„,<j)\ > nl 

Since \Br{H„, a)\ > |i?r(S„, a)\, we can replace the volume 
in &n with the volume in H„ in the last inequality. Viewed as 
a packing of Xn, the code V will then have Kendall distance 
at least to/2. ■ 

Below we consider only spheres in the space iJ„ and omit 
the reference to it from the notation Br{Hn, ■)■ 

Lemma 3.5: Let 1 = (1,1,...,!) g 7J„. Then for any 

z,y e Hn, 

2-"|i?,(z)| < |B,(1)| < \Br{y)\. 

Proof: Suppose that x ~ {xi,X2, ■ ■ ■ ,Xn) G Br{l) and 
1 ^ y = {yi,y2, ■ ■ ■ ,yn) e Hn. Consider the mapping C : 
Br{l) Br{y) where x ^ where u = {ui,U2, ■ ■ ■ ,Un) 
is given by 

jvi + {xi ~ 1) if + ix^ - 1) < n 
yn - {xi - 1) if j/j + (xj -1) > n. 

Clearly u G iJ„ and xi — 1 > \ui — yi\ for i = 1, . . . , n, 
so every point within distance r of 1 is sent to a point within 
distance r of y. Furthermore, this mapping is injective because 
if Xi,X2 are two distinct points in Br{l) then their images 
can coincide only if in some coordinates 

Vi + {xi,i - 1) = n - {x2A - !)• 

However, the left-hand side of this equality is > yi while the 
right-hand side is < yi by definition of Ui. This proves the 
right inequaUty. 



By taking the absolute values of the coordinates, any point in 
Dr{z) is sent to a point in and no more than 2" points 
have the same image under this mapping. This proves our 
claim. ■ 

niinxeif„ \Bt{Hn; x)\ ' These arguments give rise to the next proposition. 
Proposition 3.6: 

71' 71^ 
2d=T7r. :<Mn,d)< ^^ ^, ^, (11) 



where 



n\ (n-\-r—'ni—\\ 



andX„,,(i) = 

This claim is almost obvious because, by the previous lemma, 

< |C|< 



Next, 



S 

|i?,(l)| = ^Q(n,r), 



r=0 



where Q{n,r) is the number of integer solutions of the 
equation 

Xi + X2 + ■ ■ ■ + Xn ^ r, 

where < xi < n — 1, 1 < i < n. The expression for Q{n, r) 
given in the statement is well known (e.g., [13, p. 1037]). 

Expression (fTTI) gives little insight into the behavior of 
the bound. In the remainder of this section we estimate the 
asymptotic behavior of this bound and derive an estimate of 
the code capacity. 

Lemma 3.7: Suppose that r < 71^/ Inn. Then 

<Q{n,r)<i 

Proof: Let S{nJ) = E,>j (-l)'-f^n,r(i)- The lemma 
will follow if we prove that 

S'(n, 1) < and S'(n, 2) > 0. (12) 

Under the assumption on r we have 

r — ni — n + j 



fn - 

























/'r+n-n(i+l)-l\ „_i 
V r-n(i+l) / 



n 



r — ni ^ j 



/r-^n—ni— 1\ 
\ r—ni } J — 1 

i=i 
<il- 



r — n{i — 1) — 1 > 



5 



Tl-1 

< n 
V2 



< 



n 



Thus for i > 1 



Kn,r{i + 1) ^ n - « V2 

< - — ■ < 1. 



Kn,r{i) 



i + 1 n 



Therefore — iir„ ,.(2m— l)+iir„ ,,(2m) < for all m. Since the 
sum S{n, 1) starts with a negative term and the sum S{n, 2) 
with a positive one, the required inequalities in (fT2] i follow. ■ 
From the foregoing arguments we now have the following 
expUcit bounds on A{n, d) : 



n< 



2" 



^ < A{n,d) < 



J2r=0 (( 



n+1 — 1\ _ j^/r-l 



„)) 



(13) 



Here the right part is obvious and for the left inequality 
we used ( fTTT l. Lemma 13.71 and the identity ^ ^ ("^O = 

Now we are ready to complete the proof of Theorem 13.11 
Assume that d = 8(n^+'^) for some < e < 1. The two 
boundary cases of & were established in the previous section. 
Let us prove the middle equality. From ( fT3]) . 



To estimate the denominator, write 



n + t- 1 
n - 1 



t - 1 
t — n 



,t ~ n) \ t 



(because of hi(l + n/t) > (n/t) — [n/tY .) So starting with 
some n we can estimate the denominator below by V2("^l7^) ■ 
Therefore, 



2n" 



A{n,d)<—^^ 



< 



212" (n- 1) 



Next 

In A{n,Q{n^+')) 



n\nn 

On the other hand, using 



<2-(l + e) + o(l) = l-e + o(l). 



2d - 1 



< 



(n + 2d)e 



2d-\ J \ n 
and n\ > (n/3)", we obtain from (fT3T l 

n" 



)n 
< (2e2)"e(n'") 



(12e2)"e(n'")' 



Taking the logarithms and the limit, we find that "^(cZ) > 1 — e. 
This completes the proof of Theorem 13.11 



D. Bounds from embedding in Hamming space 

Since the embedding of X„ into the Hamming space Hn of 
dimension N = (2) is isometric, the known results for codes 
correcting Hamming errors can be used to derive estimates 
and constructions for codes in the Kendall space. In particular, 
the known bounds on codes in the Hamming space can be 
rewritten with respect to the space X„. For instance, the 
Plotkin bound implies that 

A{n,d) < 2d/{2d-N) 

and thus any code C C X„ with distance greater than the 
average (i.e., ^N) satisfies \C\ = 0{N). 

Given the image of a code C C Xn in Hn it is easy to 
reconstruct the code C itself. Indeed, it is immediate to find 
the inversion vector of a permutation cr given the image of a 
in Hn, and then to recover a from its inversion vector. 

Of course, not every code in Hn will have a code in 
Xn corresponding to it. The next simple proposition shows 
that nevertheless, binary codes in Hn can be used to claim 
existence of good rank permutation codes. 

Proposition 3.8: Suppose that there exists a binary linear 
[(2),/!;] code A. Then there exists an (n, > ^^Ik , d) rank 
permutation code. 

Proof: One of the 2^^*^ cosets of A in Hn must contain 
at least n\/2^~^ vectors that map back to valid permutations. 

■ 

For example, let us assume that the value N is such that 
there exists s binary t-error-correcting binary BCH code of 
length N (if not, we can add zeros to a shorter BCH code). 
Its dimension is at least N — t\og2{N + 1). This shows the 
existence of a i-error-correcting rank permutation code of size 



On the other hand, by the sphere packing bound in Hn the 
size of a t-error-correcting code in Xn is at most M < O(^). 
Thus, using the embedding Xn ^ Hn not able to close 

a gap between the existence results and the upper bounds. In 
the next section we use a different method to construct codes 
that achieve the sphere packing bound to within a constant 
factor for any given t. 

IV. Towards optimal ^-error-correcting codes 

The representation of permutations by inversion vectors 
provides a way to construct error-correcting rank permutation 
codes. In this section we construct codes in the ti space 
of inversion vectors G„ and claim the existence of rank 
permutation codes by the inequality on the code distances (O. 

We begin with constructing codes over the integers that 
correct additive errors. Once this is accomplished, we will 
be able to claim existence of good rank permutation codes. 
Let A be some subset of Z and let A^ be the space of L- 
tuples of integers from A equipped with the ti distance (|2]l. 
A code T) C is said to correct t additive errors if for any 
two distinct code vectors x,y and any 61,62 G Z"^, both of 
weight at most t, 

X + ei ^ y + e2. 
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We assume that A and t are such that V is well defined: for 
instance, below we will take A = where s is some integer 
sufficiently large compared to t. 

If in the above definition > for all i, the code is said 
to correct t asymmetric errors. However below we need to 
consider the general case, focusing on a particular way of 
constructing codes which we proceed to describe. 

Definition 4.1: Let m > L and let hi, ... , h^, < hi < 
m,i — 1, . . . , L be a set of integers. Define the code as 
follows: 



L 

C ^ e A^^^^hiXi = mod m|. 
=1 



(14) 



This code construction was first proposed by Varshamov and 
Tenenholtz [20] for correction of one asymmetric error (it was 
rediscovered later by Constantin and Rao [7] and, in a slightly 
different context, by Golomb and Welch [11]). Generalizations 
to more that one error as well as to arbitrary finite groups were 
studied by Varshamov [19], Delsarte and Piret [8], and others; 
however, these works dealt with asymmetric errors. Below we 
extend this construction to the symmetric case. 

Proposition 4.2: The code C defined in ( fT4l i corrects t 
additive errors if and only if for all e e Z^, ||e|| < t the 
sums X^i^i ^i^i ^1"^ distinct and nonzero modulo m. 
This proposition is obvious as it amounts to saying that all the 
syndromes of error vectors of weight up to t are different and 
nonzero. 

We will need the following theorem of Bose and Chowla 
[2]. 

Theorem 4.3: (Bose and Chowla) Let g be a power of a 
prime and rn = (q'^^ — l)/{q—l). There exist q + 1 integers 
jo — - ■ ■ ,jq in Z„i such that the sums 

Jii + Ji2 + ■ ■ • + Jit (0 < ii < 12 < ■ • • < it < q) 

are all different modulo m. 

This theorem provides a way of constructing an asymmetric 
t-error-correcting code of length q. This is because for any 
error vector e with ||e|| < t < m such that > 0, the sums 
Xli=i ^iji involve at most t of the numbers ji and thus are 
all different. This theorem was previously used to construct 
codes in a different context in [12] as well as in some later 
works. 

Theorem 4.4: For 1 < i < q + 1 let 



hi = 




1 + ^-2^m for t odd 
m for t even 



where the numbers ji are given by the Bose-Chowla theorem. 
Let mt — t{t + l)m if t is odd and mt — t{t + 2)rn if t is 
even. For all e G Z'+^ such that | |e| | < t the sums J^i^i ^i^i 
are all distinct and nonzero modulo nit. 



Proof: Let t be odd and let H = {0, /ii, 
Observe that 

{t - l)m/2 < h^ < {t + l)m/2. 



}• 



(i) For any ki < k2 < . . . < kt H , the sums J2i=i 
all distinct modulo m and therefore also modulo mt. These 



sums are also nonzero modulo m except for the case when all 
the fci's are 0. 

(ii) Moreover, for any ki < k2 ■ . ■ < k2t ^ H, the sum 

2t 

1=1 

and is therefore nonzero modulo mt. 

(iii) Finally, for any < fci < fc2 < . . . < fc2t G -ff and any 
r <t, 

/ ki < r — - — m < [2t — r) — - — m 



i=2t-r+l 



2t-r 
i=l 



Let us now suppose now that there exist ei, 62 G Z''+^ both 
of weight at most t such that 



9+1 

either euhi 

i=l 



9+1 



9+1 



or ^ euhi = ^ e2ihi 



However assuming this contradicts properties (i)-(iii) above. 
The claim for t even is proved in an analogous way. ■ 
Together with Proposition 14.21 this theorem implies the 
existence of a i-error-correcting code C of length n — 1 over 
the alphabet Z„ that corrects t additive errors. Recall that our 
goal is to construct a f-error-correcting code over the set of 
inversion vectors G„ which is a subset of Z"^^. Since C is 
a group code with respect to addition modulo mt, its cosets 
in ZJ^^^ partition of this space into disjoint equal parts. At 
least one such coset contains AI > n\/mt vectors from Gn. 
Invoking ^ we now establish the main result of this section. 

Theorem 4.5: Let m = ((n — 2)*+^ — l)/(ri — 3), where 
n — 2 is a power of a prime. There exists a i-error-correcting 
rank permutation code in 6„ whose size satisfies 



M > 



n\/{t{t + l)m) odd) 
n\/{t{t + 2)m) (t even). 



This theorem establishes the existence of codes whose size is 
of the same order 0(n!/n*) as given by the sphere packing 
bound of the previous section. The loss of a constant multiple 
of the optimal code size is due to the fact that we construct 
codes over the integer alphabet instead of a more restricted 
alphabet Z2 x . . . x Z„. 

As a final remark, note that the construction is explicit 
except for the last step where we claim existence of a large- 
size code in some coset of the code C. 

Acknowledgment. The authors are grateful to Gregory Ka- 
batiansky for a useful discussion of this work. 
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